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ABSTRACT 

Context. To further understand the processes involved in the formation of massive stars, we have undertaken a study of the gas 
dynamics suiTounding three massive star forming regions. By observing the large scale structures at high resolution, we are able to 
determine properties such as driving source, and spatially resolve the bulk dynamical properties of the gas such as infall and outflow. 
Aims. With high resolution observations, we are able to determine which of the cores in a cluster forming massive stars is responsible 
for the large scale structures. 

Methods. We present CO observations of three massive star forming regions with known Hll regions and show how the CO traces 
both infall and outflow. By combining data taken in two SMA configurations with JCMT observations, we are able to see large scale 
structures at high resolution. 

Results. We find large (0.26-0.40 pc), massive (2-3 Mo) and energetic (13-17 xlO'*'* erg) outflows emanating from the edges of two 
Hll regions suggesting they are being powered by the protostar(s) within. We find infall signatures in two of our sources with mass 
infall rates of order 10"* Mq yr"'. 

Conclusions. We suggest that star formation is ongoing in these sources despite the presence of Hll regions. We further conclude that 
the source(s) within a single Hll region are responsible for the observed large scale structures; that these large structures are not the 
net effect of multiple outflows from multiple Hll regions and hot cores. 
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1. Introduction 



Stars: Protostars - ISM: Hll Regions - ISM: Jets and Outflows - ISM: Kinematics and Dynamics 



The formation mechanism responsible for low mass stars like 
the Sun is, in broad terms, well constrained, from pre-stellar 
cores (Ward-Thompson et al. 2002) to pre-main-sequence stars 
(i.e. iShu et al...l987;,Andre et al...2000) . At approximately 8 Mq 
however, the situation becomes more complicated. The form- 
ing star begins to bum hydrogen hot and fast enough for ul- 
traviolet radiation to push outwards and ionize the surround- 
ing gas (although this may be shif ted to higher starting raasses 
by puffing up of th e protostar, see iHosokawa & Omukaill2009t 
lYorke & Bodenhe imer 2008|. With these extreme outward pres- 
sures, how then can accretion onto the forming star continue be- 
yond the formation of an HII region? Do es it halt before the HII 
region expands to observable sizes (i.e. lGarav&Liza"nolll999t 
iKlaassen et al] [2006)? Does it occur in self shielding toroids 
around the mass ive protostar (i.e. Keto 2007; Peters et al. 2010b; 
iLiu et al. "2010b)? Or, does a molecular accretion flow transition 
into an ioni zed one at the boundary of the HII region (i.e. lKetol 
I2OO2I l2003h ? If accretion (either molecular or ionized) contin- 
ues beyond the formation of an HII region, there should be ob- 
servational evidence for this. It could come in the form of ion- 
ized/molecular infall, or large scale outflows that are still being 
powered. 

Outflows act as an important release mechanism for the 
buildup of angular momentum in star forming regions on all 
mass scales dRicher et al.lll999tlArce et al.ll2007l) . They can also 



be some of the largest scale phenomena associated with star 
formation, and are therefore easier to observe at the large dis- 
tances to massive star forming regions. The outflows observed 
in nearby low mass star forming regions can be quite collimated, 
while coUimation factors in higher mass star forming regions 
do not appear to be as high. This last point has led some au- 
thors to suggest that as inter mediate and high m ass stars evolve 
up the main sequence (i.e. iKeto & Wood 12006). the outflows 
lose coUimation because they are older and the outflow has had 
more of a chance to wide n into the surrounding medium (i.e. 
iBeuther & Shepherdl2005l) . Indeed, there are no observations to 
date of highly collimated outflows for sources with luminosi- 
ties g reater than 10 ^ L© (masses greater than 07 main sequence 
stars: iMcKee & Ostrik er 2007). It should be noted that the high- 
est mass stars are intrinsically further away, but that the loss of 
coUimation with increasing source luminosity appears not to be 
due to beam dilution (i.e. Beuther& Shepherd 2005). Through 
characterizing the kinematics of the gas surrounding the forma- 
tion sites of massive stars, we can begin to constrain the forma- 
tion mechanisms which can produce these large scale phenom- 
ena (i.e. how much gas is accelerated into the outflow). Once 
gas kinematics are determined for a number of sources, we can 
begin to constrain models of massive star formation (i.e. give re- 
alistic parameter space boundaries for phenomena on observable 
scales). 

The large average distances to massive star forming regions 
and the clustered and embedded nature of massive star forma- 
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tion requires that these regions are observed with interferometers 
at sub-milhmeter and longer wavelengths in order to spatially 
resolve the kinematics of the gas in which massive protostars 
are embedded. With interferometers, we can resolve the spatial 
structures and kinematics of a number of molecular species and 
often determine whether the spectral line signatures are coming 
from the massive protostar/HII region or neighboring hot cores. 
In general, interferometers filter out emission from large scale 
structures and can be quite useful for ob serving small scale struc- 
tures within larger scale emission (i.e. lEisner et"aDl2008l) . With 
the Submillimeter ArrajQ (SMA) in its extended configuration, 
at 230 GHz, when structures are smaller than ~ 10", there is no 
need to include large scale structure information from a more 
compact configuration or single dish telescope since most of the 
flux is recovered by the interferometer. 

In the general interstellar medium, '^CO is very abundant 
and can be excited into some of its lowest rotational energy 
states (i.e. J=l or i-2) at the faMy low temperatures that can be 
reached in molecular clouds. Bulk motions of the gas, such as 
the material entrained into an outflow from a forming star, can 
thus be traced on large scales by CO emission. With our com- 
bined interferometer and single dish datasets, we can observe 
these large scale structures at high resolution. 

Here, we present Imm observations of '^CO and '^CO to- 
wards three massive star forming regions; GIO. 6-0.4, G28.20- 
0.05, and NGC 7538 IRSl (hereafter G10.6, G28.2, and NGC 
7538). Our observations were obtained with the SMA in its ex- 
tended configuration in late 2005, with compact configuration 
observations taken in 2007 and 2008. Single dish '^CO obser- 
vations were obtained at the James Clerk Maxwell Telescope 
(JCMXfl in 2006 (project M06BC09). In Section |2] we present 
our observations and image combination method, while we 
present the results of our kinematic studies in Section |3] In 
Sections|4]and|5]we discuss our results and conclude. 

2. Observations 

Because CO is one of the most abundant molecules in the inter- 
stellar medium (ISM), its emission is generally quite extended 
and much of the flux can be filtered out by an interferometer. 
Thus, we observed it in two SMA configurations (extended and 
compact) to observe structures on size scales ranging from ~ I" 
to ~ 25". CO emission may be more extended than this, and to 
ensure that the large scale outflows are properly characterized, 
we obtained maps of these regions at the JCMT. 

The SMA extended configuration data were obtained in 
September 2005 with CO (J=2-l, 230.538 GHz) in the upper 
sideband at a uniform spectral resolution of ~ 0.53 km s"'. The 
extended configuration data was reduced using th e SMA version 
of MIRIAD (as was done in Kla assen et alJl2009l) . 

Our SMA compact configuration data was obtained in 
September 2007, June 2008, and July 2007 (respectively for the 
source order in Table [T}. The spectral resolution was the same 
for G10.6 and G28.2 (0.53 km s '), but only half for NGC 7538 
(1.06 km s '). Just like the extended configuration data, band- 
pass calibration was done using 3c454.5, while gain calibration 

' The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 

^ The James Clerk Maxwell Telescope is operated by The Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. 



was done using the sources listed in Table |2] The compact con- 
figuration data were reduced using the SMA native MIR IDL 
package, exported to MIRIAD and self-calibrated before being 
combined with the extended configuration data. For NGC 7538, 
a two point mosaic was created in the compact configuration in- 
stead of a single pointing as was done for the other sources. Their 
pointing centers are RA (2000) = 23M3™43\75, DEC (2000) = 
-i-61°28'21".49 and RA (2000) = 23*'13'M6^75, DEC (2000) = 
61°27'59".49. 

JCMT observations were obtained for '^CO in October 2006 
using RxA. 7x7 pixel raster maps were made at 10" spacings 
at a spectral resolution of 0.63 km s ' over the full 1 GHz 
JCMT bandwidth. Each map took one hour to complete and the 
rms noise levels for each map are listed in Table |2]in units of 
Jy beam ' regridded to the resolution of the combined SMA 
datasets. At the native resolution of the JCMT observations the 
rms sensitivity hmits were 0.09 K in the T^ scale at a spec- 
tral resolution of 0.63 km s The main beam efliciency for 
these observations was 0.69. For each source, multiple JCMT 
raster maps were coadded and linear baselines removed from 
each pixel. Data reduction was completed using the JCMT star- 
link software and the datacubes were exported to MIRIAD for 
combination with the SMA observations. 

2.1. Combining Datasets 

During the formation of stars of all masses, the angular momen- 
tu m releasin g outflo ws can become quite extended. In the study 
of ' Wu et al.l (l2004l) . only two sources (one classified as high 
mass, and one classified as low mass) have angular sizes of less 
than 6" (see their Figure 3), and at most, six of the 397 sources in 
their study would have observable in the < 10" scales to which 
our highest resolution SMA observations are sensitive. Thus, to 
properly understand outflow properties at high resolution, we 
need to combine high and low resolution datasets to ensure we 
are probing the gas at all relevant spatial scales. Comparing the 
integrated flux levels in our CO observations from the SMA to 
those from the JCMT suggest that our high resolution (~ 1") 
extended configuration observations recover ~ 5% of the total 
flux and our lower resolution (~ 4") compact configuration data 
recover ~ 40% of the flux. When combined, the SMA datasets 
appear to recover approximately 50% of the single dish flux. 

Linear baselines were removed from the JCMT raster maps 
as part of the data reduction process. Thus, the continuum emis- 
sion must also be removed from the SMA datasets before combi- 
nation with the single dish data. This was done in the UV plane. 
The two SMA datasets (for a given source) were then simultane- 
ously brought into the image plane using the INVERT command 
and specifying a spectral resolution of 0.63 km s"' to match the 
spectral resolution of the JCMT observations. For NGC 7538, 
the SMA compact configuration data had a spectral resolution 
of 1.06 km s so, when inverting, this lower spectral resolution 
was used instead. 

Since combining the JCMT data with that from the SMA 
is a form of cleaning, both a JCMT image and beam map are 
required. The JCMT beam map used for deconvolution was cre- 
ated using the MIRIAD command IMGEN assuming a gaussian 
source of intensity 1 and a radius of 20.8" (the resolution of 
the JCMT at 230.5 GHz). The image cube was then regridded 
to the resolution of the jointly inverted interferometer image. 
Once the single dish and interferometer maps were on the same 
flux scale and at the same resolution, they were combined us- 
ing a maximum entropy method (using the MIRIAD command 
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Table 1. Observed Sample of Massive Star Forming Regions. 



Source 


pointing center (J2000) 


Vlsr 




Distance 


RA DEC 


(km s"') 


ref 


(kpc) ref 


GlO.6-0.4 


18:10:27.77 -19:56:04.5 


-3 


1 


6.0 4 


G28.20-0.05 


18:42:58.17 -04:13:57.0 


99 


2 


5.7 2 


NGC7538 IRSl 


23:13:45.30 +61:28:10.0 


-59 


3 


2.7 5 



References. (D lPurcell etai] ( 120061) : (2) lPestalozzi et all ( 120051) : (3) IZheng et all ( 120011) : (4) lKeto et al] 120081) : (5) lMoscadelli et all ( 120091) 
Table 2. SMA and JCMT Observation Characteristics 



Source Gain Calibrators Res." RMS noise (Jy/beam) Synthesized Beam'^ 

extended compact (kms"') SMA'' JCMT" Comb? Majo? Minor" PA" 

Gloe 1833-210 1733-130 (163 024 09 024 TaP UlV' 54A^ 

G28.2 1944-201 1733-130 063 012 1.2 014 1.88" 1.07" 77.9° 

NGC7538 2202+422 0102+584'' 1.06 Oil 2.0 013 2.35" 1.61" 71.6° 



Notes. 

Spectral resolution of the final data products. Synthesized beam of the combined SMA and JCMT dataset. Jointly inverted SMA 
extended and compact configuration data. JCMT observations regridded to the spatial resolution of the jointly inverted SMA observations. 
'''' RMS noise value for the combination of the JCMT data and the jointly inverted SMA data using the miriad task MOSMEM at the spectral 
resolution of the final data product. BLLac was also used as a gain calibrator, and 3c273 and Uranus were both used as bandpass calibrators in 

addition to 3c454. 



MOSMEM) and a clean image was restored using the interfer- 
ometer map as the dirty image and the results from MOSMEM 
for the clean components. 

For each source, we also tried combining the single dish and 
interferometer maps by feathering (using the MIRIAD command 
IMMERGE) which did not recover all of the single dish flux. 
We also tried inverting the JCMT observations to the UV plane 
and jointly inverting all three datasets; however this method only 
produced reliable results for two out of the three sources, and for 
consistency, we present our results of the MOSMEM combined 
datasets for all three sources. To test the robustness of our image 
combination, we determined outflow properties from the single 
dish only data, the single dish and lowest resolution interferom- 
eter map, or all three datasets combined. Using the same spa- 
tial and velocity limits, the properties are the same in all three 
instances. The values are lower however when we only com- 
bine the two interferometer maps because we are filtering out 
the large scale structures consistent with the flux filtering de- 
scribed earlier in this section. We have also compar ed ou r JCMT 
obsei-vations of NGC 7538 to those of iPavis et all (Il998l) . Using 
their velocity integration limits (-94 to -64 km s~^ and -50 to -20 
km s"'), we compared the fluxes in our map to theirs. Our map 
is much smaller than theirs, and they show CO contours well be- 
yond our map boundary. Their emission covers a much greater 
area than we have sampled with the SMA, and thus we can re- 
ally only compare peak fluxes between the two maps. We find 
that our fluxes are somewhat lower than theirs, however both the 
red and blue peaks are within 15% of theirs. 



3. Results 

Here, we present our observational results. Of the three observed 
sources, we find indications of large scale molecular outflows in 
two sources (NGC 7538 and G28.2), as well as molecular infall 
signatures in two sources (NGC 7538 and G10.6). We note that 
while there is no large scale molecular outflow signature coming 
from G10.6, an ionized outflow has been previously detected in 
this region (iKeto & Woodll2006h and despite our non-detection 



of infall i n G28.2, an infall sig nature has been seen in the dense 
gas (NH^. ISollins et alll2005bl) surrounding this region. 

3.1. Outflow of molecular gas from the HII regions 

While outflows are generally characterized as one of the larger 
phenomena associated with star formation, it is only by observ- 
ing them at high resolution that we can constrain which source is 
producing the observed large scale outflow. In each of the cases 
presented below, when detected, the outflow appears to emanate 
from a single HII region. The molecular outflow originates at the 
boundary of the HII region, suggesting that perhaps the outflow 
generation is occurring within the HII region. 

The three datasets for CO were combined as described in 
Section im and from these datacubes we created zeroth moment 
(integrated intensity) and first moment (intensity weighted ve- 
locity) maps of the red and blue shifted emission independently 
for each source to determine the outflow dynamics. These maps 
were created using both velocity range and noise clipping masks. 

The high abundance of CO does make it an excellent tracer 
of high velocity and bulk gas motions; however, it also causes 
the line to become optically thick quite quickly. For our calcula- 
tions we have only integrated over the high velocity gas regions. 
If we were to integrate over the lower velocity gas, we begin de- 
tecting ambient and envelope gas. Our inner velocity limits were 
chosen though trial and error. We started with our low velocity 
limits quite close to the Vlsr, and continued moving that limit 
to higher blue and red shifted velocities until the maps showed 
primarily outflow (i.e collimated) structures, with little ambient 
contamination. Our inner velocity limits are at relatively high 
velocities, which means we are only integrating over the 'outer 
wings' in the terms of (Margulis & Lada 1985). This is a limi- 
tation of our observations, and thus we are underestimating the 
outflow properties. We have assumed optically thin emission in 
the wings. 

In all of our calculations below, we assume an ambient tem- 
perature of lOOK. This value was chosen to represent the bulk 
of the outflowing gas since it is colder than the hot core gas 
surrounding each of these HII regions ( > 300 K as traced by 
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(a) G28.2-0.04 
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23''13'»48* 47* 46* 45" 44* 
Right Ascension (J2000) 

(b) NGC 7538 IRS 1 



Fig. 1. Large-scale Molecular Outflows: In both figures, the blue and red contours show the blue and red shifted CO emission 
integrated over the velocity ranges given in Table |2] starting at 5cr (3cr), and increasing in intervals of Icr. The green contours 
represents the continuum emission, starting at 50% of the peak emission. In the right panel, the numbers indicate the locations of 
IRS 1,2, and 3. 



CH3CN, iKlaassen etal1[2009h and warmer than the tempera- 
tures found in cold dense cores (< 20 K i.e. iGarav et al.ll2004t 
iPeretto et al.l l201d). We note that changing the temperature by a 
factor of two changes our estimates by less than a factor of two. 

Integrated intensity maps were created for the sources with 
outflows and are presented in Figure [1] while the velocity ranges 
used to create the maps (centered on the local standard of rest 
velocity (Vlsr) given in Table [1) are presented in Table |3] To 
create the first moment maps and determine the average red and 
blue shifted velocities, we shifted the velocity axis of each dat- 
acube so that the Vlsr is centered at km s"', and set range 
masks over the same velocity ranges used for the zeroth moment 
maps and clipped the images. 

Using the spatially averaged integrated intensity in both the 
red and blue outflow lobes, we calculated the average H2 col- 
umn density of the emitting region by using a constant of pro- 
portionality between the column density and integrated inten- 
sity (N= 1.3 X 10'^ jTdv) which assumes a temperature of 
100 K. This average column density, when multiplied by the 
emitting area and mass of H2 gives the mass in each outflow 
lobe (M = NA) when the abu ndance of CO with respect to 
H2 is taken into account (10~^. iHerbst & vaiT Pishoeck 2009). 
Since a first moment map is in essence a measure of veloc- 
ity, the zeroth and first moment maps were multiplied together 
to determine the momentum (P - Mv) in each outflow lobe. 
Similarly, we also determined the kinetic energy in each lobe 
(E - 1/2 * Mv^ = P^/2M). These mass, momentum and energy 
values are presented in Table |3] The quoted uncertainties come 
from the uncertainties in the integrated intensities and velocities, 
and are propagated through the equations described above. 

We estimated a kinematic age for G28.2 (tdyn - 2.4 x lO'* yr) 
by dividing the size of the outflow by the velocity of the gas. For 

NGC 7538, we took a value from the literature (tdvn = 4.2 x 10"* 

I " * 

vr. Davis et al. 1998) since larger scale observations of this re- 
gion show outflow motions well beyond the edges of our pri- 
mary beam. Using these kinematic ages, and the derived out- 
flow masses and energies, we determined average mass loss rates 



Table 3. Outflow Results from Kinematic Analysis 



G28.2-0.04 
(fdyn = 2.4 X 10* yr) 
Blue Red 
(Av=-10,-3) (Av=3,10) 


NGC 7538 IRS 1 
(tay„ = 4.2 X 10* yr) 
Blue Red 
(Av=-20,-10) (Av= 10,20) 


M" 


0.89 


1.60 


0.43 


0.98 


P* 


8.3 


8.7 


4.6 


10.7 




7.6 


4.9 


4.9 


11.8 


L'' 


0.30 


0.20 


0.10 


0.23 




3.6 


6.4 


1.0 


2.3 



Notes. A calibration uncertainty of 20% is not propagated through these 
equations. The velocity ranges listed above each column (Av) represent 
the integration limits with respect to the Vlsr stated in Table[T] 

Mass in outflow lobe in units of Mq Momentum in outflow lobe 
in units of M© km s"' Energy in outflow lobe in units of 10** erg 
Mechanical Luminosity in outflow lobe in units of Lq Mass 
outflow rate in units of 10"' M© yr"' 



(M = Mjt) and outflow mechanical luminosities (L = £/?)• 
These values can also be found in Table [3] For both outflow 
sources, the mass outflow rates are greater than the values gen- 
erally seen in lower mass star for ming regions (w hich can be 
in the range 10-8 - 10-5 Mo yr-\IShepherd"2003'). Assuming 
the amount of mass in the outflow scales proportionally to the 
amount of mass infalling onto the star (Beu ther et al.t 2002). then 
we would expect larger mass outflow rates than for lower mass 
stars since the masses of the outflows are larger. 

For each source, the orientation of the outflows are gener- 
ally consistent with being perpendicular to the rotation directions 
of the hot cores surrounding the HII regions (as traced by OCS 
and SO2 in lKIaassen et alJ[2009l) . and t he outflow kinematics are 
consistent with other outflows listed in IWu et al.1 (12004 and to- 
wards the values they show for high mass star forming regions 
(those shown with open squares in their figures). We note that 
one of our sources overlaps with their study (NGC 7538), and 
that our derived outflow properties are consistently lower than 
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their values. As discussed in Section |2] our map does not cover 
the entire flow in NGC 7538, and we are thus missing much of 
the extended flux. Since mass scales directly with integrated in- 
tensity under the assumption of optically thin gas, we expect our 
outflow masses to be lower than p reviously derived values (i.e. 
iDavis et iDIIMl lOiu et alj|20Tll who used larger JCMT maps 
to determine the outflow properties). It is likely that this is the 
case for G28.2 as well. 

The lower derived outflow properties can also be understood 
in terms of the high op tical depths of CO at low velocities. 
iMar^uUs & Ladal (Il985l) discusses how the bulk of the outflow 
mass and dynamics can be found in the lower vel ocity wing 
emission, from what they called the 'inner wing'. IChoi et al] 
(Il993h described this same region as containing the 'high veloc- 
ity' gas. To use the terminology from these papers, we are only 
tracing the gas in the 'outer wing' or 'extremely high velocity' 
regions (respectively). As stated previously, '^CO becomes opti- 
cally thick quickly in massive star forming regions, and we can- 
not trace the bulk of the gas that is optically thick in this tracer 
The high opacity of the line is advantageous for seeing the ex- 
tremely high velocity gas, but means that we are missing most 
of the mass in the outflow. We cannot comment on the mass of 
the gas at lower velocities than the limits quoted in Table [3] but 
must assume that it is quite large if there is still more than a solar 
mass of material in the outer wings of these outflows. 



3.2. Infall of Molecular Gas onto the Hi! regions 

Bulk infall motions can manifest as inverse P-Cygni pro- 
files with reds hifted absorption and blueshifted emission (i.e. 
iDi Francesco et al. 200 1) along the line of sight to the infalling 
region against a bright background. From this profile, a rough 
estimate of the infall velocities and masses can be determined. 
[Beltran et al. (2011) found inverse P-Cygni profiles in '^'CO to- 
wards two of the three sources they studied; the third, G10.6, 
had i^CO emission on scales much larger than they were able 
to probe with their data. For the sources where they do detect 
infall, they find infall velocities of 3.8 and 5 km s ' for G19.61- 
0.23 and G29.96-0.02, respectively. Similarly, inverse P-Cygni 
profiles hav e also been observe d towards G24. 78 -1-0.08 ( 2 km 
s ' in NH. 'Beltran et al.l l2006h . G31.41-h0.31 ( 3.1km s ' in 
C^'^S Girart et al. 2009|rand in W51 IRS2 ( 4 km s ' in CN 
IZapata et al.ll2008i) . In order to determine whether this infall sig- 
nature is present towards our sample of sources, we used our 
combined JCMT and SMA CO data. We note that the infall sig- 
natures, when present, appear virtually the same whether or not 
the JCMT data is included since these signatures occur on small 
scales (a few arcseconds). We find infall signatures in the CO 
from G10.6 and NGC 7538, but not towards G28.2. Below, we 
describe how these signatures were characterized to determine 
source properties. 

We found the pixel with the largest CO absorption feature 
along the line of sight to the HII region, and created an integrated 
intensity map based on the width of the absorption feature. We 
determined the -3cr contour of the feature in the moment map, 
and created an integrated spectrum based on the overlap between 
that contour, and the surface of the HII region. For G10.6, the 
smaller region was the absorption feature, for NGC 7538 it was 
the unresolved HII region. 

U sing the two layer infall model of iDi Francesco et all 
(l20()ih we determined infall velocities from the inverse P-Cygni 
profiles in CO for G10.6 and NGC 7538. The inputs to these 
two layer models consist of Planck functions characterized by 
their different temperatures. These temperatures correspond to 



Table 4. Infall Model Parameters from Inverse P-Cygni Profiles 



Source Name 


T, 


Tc 




Vlsr 


M,„ 




(K) 


(K) 


(kms-') 


(km s-') 


(Mo yr-i) 


G10.6 


250 


170 


2.5 


1 


2x10"'* 


NGC 7538 


385 


400 


2.3 


-59 


2x10-'' 



the cold layer of gas in front of the emission region (Tf=7.4 K, 
the excitation temperature of CO J =2-1), to the rear of the emis- 
sion region (T^) and the temperature of the continuum source 
(Tf). The infall velocities are presented in Table |4] With these 
infall velocities, known source distances and integration areas, 
we also determined the mass infall rates using Equation 3 from 
iKIaassen & Wilson ( 2007) : 



(1) 



where is the ambient density, set here to 10* cm"-', fi 
is the mean molecular weight (2.35), ;«// is the mass of a hy- 
drogen atom, r is the radius of the emitting area, and v/„ is the 
infall velocity. We suggest that these inverse P-cygni profiles are 
due to local infall and not foreground absorption because these 
same structures have been seen in other molecules towards these 



sources as well. S ee for instancelQ iu et al,l(2011 ) for NGC 7 538 
and for G10.6 see lHo & Haschick ( 1986); Keto etaP ( [198^) 

For both sources, the mass infall rate was 2x10"'* Mq yr"'. 
For G10.6, the rest velocity required to align the model with the 
observations is shifted with respect to the local standard of rest 
velocity by -i- 4 km s '. The profiles plotted in Figure |2] only 
reflect the CO along the line of sight to the HII regions, and not 
the areas shown in Figure [T] 

In both sources, inverse P-Cygni pro files can also be seen in 
CH3CN (J=12-l 1) originafly presented in lKlaassen et al.l(l2009l) . 
The derived temperatures (323 and 310 K respectively) and spa- 
tial structures of the CH3CN regions in G10.6 and NGC 7538 
as described in lKlaassen et al.' suggest that this species is tracing 
the hot cores surrounding these two HII regions. The profile is 
strongest in the K=3 transition for both sources, and these in- 
fall signatures are consistent with those seen in CO; the 4 km 
s ' velocity shift in the G10.6 profile included. That the CH3CN 
also shows an infall signature suggests that material is still ac- 
tively accreting onto the two HII regions. The CH3 CN spectra 
shown in the top and middle panels of Figure 7 of Kl aassen et al.l 
('2009') are integrated over the entire CH3CN emitting region for 
G10.6 and NGC 7538 respectively, and not just the area show- 
ing an inverse P-Cygni profile, yet the signature is still visible 
in the integrated spectrum from NGC 7538 (see Oiu et al. 201 ll). 
The simulations of high mass star forming regions of Wa ng et al.l 
(2010) which include both outflow feedback and magnetic fields 
find similar mass infall rates, and, as discussed later, this molec- 
ular mass infall rate in G10.6 may be comparable to the ionized 
mass infall rate as traced by H66q;. 

3.3. Results for Individual Sources 
3.3.1. NGC 7538 IRS1 

As shown in the right panel of Figure [T] both the red and blue 
outflow lobes appear emanate from the central HII region (IRS 
1), suggesting that the outflow may be precessing or bending 
around other structures (such as IRS 2). On much larger scales 
than observable here, the outflow appears to be dominantly in 
the North-South direction (S. Corder, priv. comm). The outflow 
direction we find is consistent with the previous observations of 
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(a) GlO.6-0.4 (b) NGC 7538 IRS 1 

Fig. 2. CO (J=2-l) spectra integrated over the infall regions overlapping with the HII regions for G10.6 (left) and NGC 7538 (right). 
The grey lines shows the 2 layer infall model characterized by the parameters given in Table|4] 



iDavis et a which were on similar size scales to those 

observed here, but at much lower resolution. 

The energetics from this source are comparable to those in 
G28.2, despite having a lower bolometric luminosity than G28.2. 
This is interesting because it is much closer than either G28.2 
or G10.6. It is unclear how distance affects our results since the 
outflow from NGC 7538 would still be resolved at any of the dis- 
tances to the sources studied here. The outflow from this source 
is the most extended in this study, with each lo be spanning more 
than .4 pc (at an assumed distance of 2.7 kpc iMoscadelli et all 
I2009h . 

Comparing our CO map to Spitzer IRAC images of the re- 
gion, we see that the blue outflow lobe curves around a 3.6 and 
5 fim feature coming from IRS 2 which is not very prominent in 
the 4.5 jum band. This could mean that the CO is tracing the edge 
of a region that is deficient in molecular gas, and dominated by 
PAH emission since H2 and CO are expected to emit in the 4.5 
fj.m band while PAH emission is expecte d to be more prom inent 
in the 3 and 5 jum bands (see Figure 1 of iReach et al.ll2006h . The 
portion of the Spitzer map corresponding to the red emission 
lobe is saturated in all four IRAC band s. The ioniz ed emission 
from this region (as traced by H30q' see lKeto et ai][2008l) IS un- 
resolved in these observations and does not show any systematic 
or bulk motions, so a similar analysis cannot be conducted. 

Because NGC 7538 is relati vely nearby for a high mass star 
forming region (2.65 kpc Mosc adelH et"aLll2009l) . a number of 
distinct sites of massive star formation have been identified (i.e. 
IRS 1, 2, and 3). Both IRS 2 and 3 are within our field of view, as 
indicated in Figure [T] yet it is clear that the detected large scale 
outflow is centered on IRS 1 . 

As stated earlier, the outflow properti es derived here ar e 
much smaller than thos e determine d in iDavis et al ] (Il998h . 
Judging from Figure 12 of lDavis et al.L their intensities were cal- 
culated over an area greater than 80" X 80", while our integration 
area was closer to 40" x 40" since we only had a single point- 
ing wi th the SMA in its extended configuration. iKameva et alJ 
(119891) observed CO J=l-0 in this source with flie NRO 45 m 
antenna. Their map of the gas coming from NGC 7538 IRSl 
is only slightly more extended than ours, and their derived out- 
flow properties are closer to our values than those of Davis et al.l 
The maps and derived outflow properties for this source from 



lOiu et al.1 (1201 ih are also larger than ours. This suggests that 
their larger maps are able to detect more of the flux fro m the 
outflow. That the outflow properties in lKameva et al.l ( 1989.) and 
iOiu et al.1 (l201ll) are lower than those of D avis et'aiT(ll998l) sug- 
gests t hey have al so missed part of the outflow. It should be noted 
that in lOiu et"al] (|2011), the outflow parameters are calculated 
over the multiple outflows they suggest for this region. 

The continuum brightness for this source at 230 GHz is 
~ 2.83 Jy/beam, and the absorption feature is strongest ( -2.9 
Jy/beam) at -54 km s"' suggesting that there is complete ab- 
sorption along the line of sight. We suggest there is complete 
absorption over a velocity range of a few km s"' (as shown by 
the model spectrum) with contaminating emission from the out- 
flow causing the line intensity to increase to near zero at -56 km 
s-i. 



3.3.2. G28.20-0.05 

The 5.8 yum Spitzer IRAC image of this region shows a slight 
extension in this region (centered on the HII region) towards the 
south east of the map in the left panel of Figure [1] It should be 
noted that this extension is perpendicular to the outflow shown 
here a nd along the rotation direction shown in iKlaassen et al.l 
(2009^), but is much more extended than the warm rotating gas 
seen in OCS and SO2. We also note that the radio recombination 
line (H53g. lKeto et al. 2008) data for this source does not appear 
to have an infall signature, but there does appear to be a velocity 
gradient across the source consistent with the rotation seen in 
IKlaassen et alj ( 120091) . 

As was the case for NGC 7538, the outflow moment map- 
ping appears to fit best when the LSR velocity is taken from the 
average of the values listed in the literatur e over the last ten year s 
(99 km s ') instead of the one assumed in IKlaassen et aT] ( l2009l) . 

The outflow kinematics from this source are all sl ightly 
lower than the trend shown in lLopez-Sepulcre et al.l(l2009L sug- 
gesting we are perhaps missing some of the gas at high opti- 
cal depths. The red and blue outflow lobes span 0.55 and 0.41 
pc on the sky (respectively, at an assumed distance of 5.7 kpc) 
which makes the physical scale of this outflow comparable to 
NGC 7538. 
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Curiously, we do not detect an infall signature towards this 
HII region like we do in the others, despite the infall signature 
having been seen in NH3 by Sollins et al. (2005a). This may sug- 
gest the large scale infall has stopped but that infall is still ongo- 
ing in the densest gas surrounding the HII region. 



3.3.3. G1 0.6-0.4 

An ioniz ed outflow has been seen towards G10.6 in H66a 
dKeto & Wood .2006). but this is the first high resolution search 
in this source for a molecular outflow which includes the large 
scale structures observable with a single dish telescopes. That 
there is no evidence for a Jarge scale m olecula r outflow is highly 
curious (although see iUu et al]r201(j3 1201 ll for a possible ex- 
planation). Previous studies of this region have suggested an ori- 
entation that is primarily in the plane of the sky. 

Inverse P-Cygni profiles exist in both CO and CH3CN for 
this source. As shown in Table |4] the two layer infall model 
is only able to reproduce the CO inverse P-Cygni profile if the 
rest velocity is set to 1 km s~' (not -3 km s"' as listed in Table 
[U. That these two species show this absorption feature slightly 
redshifted from the source velocity suggests that perhaps the 
molecular gas is doing something different than the ioni zed gas. 
It shou ld be noted that the NH3 observations of (Solhns et al.l 
I2005b^ for this source also show a similar velocity shift at high 
resolution for the molecular gas. The first moment map they 
present in Figure 2 shows gas between -3 and 4 km s"' suggest- 
ing an average velocity of 1 km s ' as well. IKeto et alj ( 1198^ 
also place the central velocity of their NH3 observations at 1 km 
s"' (see their Figure 1). 

There is the posibility that our infall signature is contami- 
nated by absorption from the foreground spiral arm in the ap- 
prox imate velocity ra nge of -10 km s"' to 50 km s"' (see Figure 
1 of EiE et al.ll2003h . However, that absorption feature is quite 
flat in HI and thus we suggest it only minimally contributes to 
our absorption feature. Our feature has a depth consistent with 
total absorption over a much smaller velocity range than the 60 
km s"' range of the HI absorption sugge sting t hat this is a sepa- 
rate feature (although see Appendix A of lLiu et al. 2010a). None 
of the parameters going into our infall model are greatly affected 
by this absorption as the depth of our absorption feature shows 
complete absorption. 

The molecular infall appears to be consistent with the infall 
of the ionized gas as traced by H66a (Keto & Wood 2006) . From 
modeling the infall signature in the ionized gas, Keto & WoodI 
('2006') determined a mass infall rate of 1x10 ^Mq yr From 
our two layer infall model, we determined a molecular gas infall 
veloci ty, which we input into Equation 3 of iKlaassen & Wilsonl 
(12001) and determined the molecular gas infall rate to be 
2xl0~"^ Mo yr~', the same as th at infer r ed fro m the NH3 observa- 
tions of iHo & Haschickl ( ll986h . lKetd d 1990 1) show a mass infall 
rate of 2.5x10""* at a distance of O.I pc from the HII region as 
well. 

For comparison, m olecular infall ha s also been detected to- 
wards G24.78-i-0.08 (B eltran et al.l2006h . and that HII region ap- 
pears to be undergoing a contractio n consistent with an accretio n 
flow onto the HII region as well (iGalvan-Madrid et~al ] l2008h . 
This suggests that G10.6 is not the only source with these prop- 
erties, although G2 4.78-h0.08 has been shown to be powering a 
large scale outflow dFuruva et al.ll2002l) . 



4. Discussion 

Here, we only see one outflow emanating from each massive star 
forming region. This is not always the case however; in G5.89- 
0.39, the large scale outflow (Klaassen et al. 2006 ) is in a dif- 
ferent dire ction than the smalle r scale outflows of ISoUins et al.l 
(120041) and lHunter et al.l (l2008l) who also see multiple outflows 
in various directions (see Oiu et al.ll201 ll for an alternate inter- 
pretation of NGC7538X 

Figure 5 of iLopez-Sepulcre et al.l (l2009h shows outflow 
mass, momentum and energy as a function of bolometric lumi- 
nosity for a large sample of outflows from massive star forming 
regions observed with the IRAM 30 m telescope. They find that 
each of these quantities increases as a function of source bolo- 
metric luminosity. We find that our derived values are lower than 
their trend, which we suggest is due to using '^CO instead of 
'^CO for our analysis. We are primarily tracing 'outer wing' gas 
(in the terminology of Margulis & Lada 1985), where, by defini- 
tion, the '^CO is no longer emitting. Were we able to distinguish 
between the ambient and outflow gases at lower velocities, our 
fluxes would increase, and thus our deri ved parameters would 
be much more consistent with those of lLopez-Sepulcre et al.l 
(2009). An overestimate in the flux will lead to a proportional 
overestimate in the outflow properties (see the equations given 
in SectionjTT)- 

The advantage of the higher resolution data is that we are 
able to determine which source i s powering the large sca le out- 
flow with much higher precision. [Keto & Klaassenl d2008l) at the 
same resolution as the data presented here, showed that the small 
scale ^^CO surrounding W51e2 appears at the edges of the HII 
region (see their Figure 2). The outflow in this region can be 
traced all the way down to the ionization boundary, despite the 
presence of multiple HII region s within 10" of each other (see, 
for instance IZhang & Hol[T997l) .Thus, in order to reveal which 
source is responsible for the large scale energetics, high resolu- 
tion observations are crucial. Here, we have shown that the large 
scale outflows are produced within a single HII region, and are 
not the prod uct of the energetics of the surrounding hot cores. 
Peters et al. ( 20T0ah simulated accretion and feedback in a mas- 
sive star forming region with multiple sink particles representing 
multiple massive stars. They showed that the single large scale 
outflow produced from these clusters comes from the most mas- 
sive member of the cluster. Our observations seem to be showing 
the same thing; the large scale outflow coming from the brightest 
HII region within our field of view. 

In both cases, we have been able to trace the outflow back to 
its powering source (the HII region). That the outflowing gas can 
be traced back to the edges of a single HII region suggests that 
the massive protostar(s) within that HII region are powering the 
outflow. We also suggest that something inside the HII region is 
powering the large scale structure because of the continued infall 
seen towards these sources as well. This leads us to conclude that 
the massive, large scale outflows detected in these regions are 
powered by the massive protostar(s), and are not the averaged 
gas motions from the entire cluster expelling gas to release an- 
gular momentum. Our ratio of Mout to Mi,, is approximately 0.3, 
which appears to be consistent with the ratio of the jet entrained 
outflow mass to mass accretion rate discussed in Beuth er et^ 
(2002), noting that they found a lot of scatter in their value (0.17 
for a source whose luminosity is accretion driven). 

There are infall signatures for each source, either presented 
here or in ISollins et al. (2005a) suggesting continued accretion 
beyond the formation of the HII region. In none of our sources is 
this as clear as in G10.6 where both molecular and ionized infall 
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are detected. The infalling molecular gas appears to pass through 
the ionization boundary, become ionized and continue infalling 
with similar mass infall rates. Spatially, these two infalling pop- 
ulations appear to be nested within each other, supporting this 
theory further. These observations further suggest that high mass 
star formation may be able to proceed in a manner similar to, but 
scaled up from, lower mass star formation with the added com- 
plications of HII regions and clustering. 

It is interesting that we do not detect an infall signature in 
G28.2, wher e one has pre viously been detected in the dense gas 
(in NH3, see lSoUins et al.l l2b05a). This suggests that perhaps the 
large scale collapse has finished, but infall and accretion are on- 
going in the dense gas close to the HII region which we cannot 
trace with CO or CH3CN. Oddly, this source is our youngest 
source (based on a kinematic age for the outflow) as well as the 
one with the most collimated outflow. Further investigation will 
be necessary to understand why infall was not detected in CO. 

5. Conclusions 

We combined observations of CO (J=2-l) from two SMA con- 
figurations (extended and compact) with observations of the 
large scale CO structures observed at the JCMT for three mas- 
sive star forming regions. We detected large scale molecular out- 
flow signatures from two of them as well as molecular infall 
signatures towards two sources one of which has previous de- 
tections of an ionized infall in H66a. 

Due to the high resolution of our datasets we have shown 
that the outflowing gas originates at the edges of the HII regions. 
This suggest that the outflows are driven from a single powering 
source and are therefore not the spatially averaged outflows from 
a star forming cluster. Our data also suggest that these outflows 
are continuing to be powered despite the presence of the HII 
regions. 

If these outflows are powered by accretion, which is likely 
the case, the continued presence of powered outflows is an in- 
dication that accretion is ongoing as well, despite the presence 
of an HII region. This is further supported by the detection of 
infall signatures in the molecular gas towards two of our sources 
which small scale infall has been detected in denser tracers by 
other authors for the third source. From inverse P-Cygni profiles, 
we infer infall velocities and masses much higher than for lower 
mass star forming regions, but consistent with those found in re- 
gions forming high mass stars. We also note that, for G10.6, the 
molecular and ionized infall rates for the ionized and molecular 
gas even appear to be similar. 
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